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Ab initio determination of polarizabilities and van der Waals coefficients of Li atoms 
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We report a new technique to determine the van der Waals coefficients of lithium (Li) atoms based 
on the relativistic coupled-cluster theory. These quantities are determined using the imaginary parts 
of the scalar dipole and quadrupole polarizabilities, which are evaluated using the approach that we 
have proposed in [l[. Our procedure is fully ab initio, and avoids the sum-over-the-states method. 
We present the dipole and quadrupole polarizabilities of many of the low-lying excited states of Li. 
Also, the off-diagonal dipole and quadrupole polarizabilites between some of the low-lying states of 
Li are calculated. 

PACS numbers: 31.15.Ar,31.15.Dv,31.25.Jf,32.10.Dk 
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I. INTRODUCTION 

In recent years, ultra-cold atom experiments have been 
used in the study of a variety of scattering physics, in- 
cluding the probing of different types of phase transitions 
|2J . From an experimental point of view, lithium (Li) is a 
very interesting system since its 6 Li and 7 Li isotopes cor- 
respond to fermionic and bosonic systems, respectively. 
These isotopes are used in the study of boson-boson y, |j| , 
boson-fermion 5] and fermion-fermion mixtures (J,|6|. 

For the theoretical description of these kinds of sys- 
tems, a knowledge of the interatomic potential is neces- 
sary. At a large nuclear separation R, the s-wave scat- 
tering interatomic potential is accurately represented by 
the sum of two independent contributions, the exchange 
and electrostatic potential Q- The former is related to 
the ionization energies and scattering lengths which will 
not be discussed hereafter. The electrostatic potential is 
given by @ as 



V{R) --&-& 



(1.1) 



where C§ and Cs are known as dispersion or van der 
Waals coefficients. As R — > oo, the long-range poten- 
tial V(R) is dominated by —Cq/R 6 and —Cs/R 8 , where 
the higher-order terms are sufficiently weak to be ne- 
glected. Both coefficients can be evaluated from the 
knowledge of the imaginary parts of the dynamic dipole 
and quadrupole polarizabilities [3, [l0(. Several groups 
have evaluated these quantities because of their neces- 
sity in the simulation, prediction, and interpretation of 
experiments on optical lattices in cold-atom collisions, 



photo-association, and fluorescence spectroscopy 111 

Since the classic work of Dalgarno and Lewis 13j ] . 
different procedures have been followed to deter- 
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mine polarizabilities. An often- used method is the 
sum-over-intermediate-states approach, which employs 
dipole/quadrupole matrix elements and excitation en- 
ergies of important states [Tj, [ija, [lg, [l?], [T$|. Sum- 
over-the-states methods are, however, limited in their 
accuracy because of the restrictions in the inclusion of 
higher states, which are difficult to generate. Coupled- 
cluster based linear response theory [lj, UM seems to be 
a promising method to study both static and dynamic 
polarizabilities, while avoiding this limitation of the sum- 
over-the-states approach. This method is well applicable 
to closed-shell systems. For relativistic open-shell sys- 
tems and adopting atomic symmetry properties, however, 
it is not an easy formalism. Therefore, the sum-over-the- 
states approach using dipole/quadrupole matrix elements 
or oscillator strengths is often used in open-shell atomic 
systems [THUS- 

In this work, we present a novel approach to deter- 
mine the van der Waals coefficients for lithium using a 
method which employs fully atomic symmetry proper- 
ties in the framework of the relativistic coupled-cluster 
(RCC) approach. The approach is ab initio and avoids 
the limitations of the sum-over-the-states methods. It 
has recently been employed to determine static polariz- 
abilities in closed-shell and one-valence open-shell sys- 
tems p], Ojl, E(|. We also present the static dipole and 
quadrupole polarizabilities of many of the excited states 
in Li. These could be useful in the calculation of the 
dispersion coefficients of the excited states and in deter- 
mining Stark shifts. So far, only a few studies have been 
carried out on t he p olarizabilities of the Li excited states 
EH HI [H El El [1|. Most of these studies use non- 
relativistic theories, and we will compare those results 
to our relativistic calculations to assess the relevance of 
relativistic effects. We also present the scalar polarizabil- 
ities among two different states, which are of interest for 
several types of studies [27| . 

The outline of the rest of the paper is as follows. We 
start by presenting the theory for polarizabilities and van 



der Waals coefficients in Sec. II. Next, we discuss our 
method of calculation in Sec. Ill, and in Sec. IV we 
present and discuss our results. 



II. THEORY 

In this section we give the definitions of the static and 
dynamic polarizabilities and the van der Waals coeffi- 
cients. 



A. Polarizability 

The static dipole polarizability ai(J 1 ,,mj„) of a valence 
(v) state \^ v ) of a single valence system is given by |28l . 



a>i(J v ,mj v 



(^|£>|^)(^ 7 |£>|W„) 
Ei-Ev 



E 

0/, x , 3m l ~ J«(J« + 1) 2 

ai(Jv) H 77777 7T ai(Jv), 



(2.1) 



where the scalar polarizability ai(J v ) is given by 

and the tensor polarizability af (J„) by 

10J o (2J u - 1) 1 1/2 



af(J„) 



EH*-{* i *} 



3(J„ + 1)(2J„ + 1)(2J„ + 3) 

W JJ« 1 J/] |q||J>||J/)| 2 



(2.3) 

Here £> is the dipole operator, J v and mj v are the an- 
gular momentum quantum numbers of \^f v }- \^i) rep- 
resents allowed intermediate states with respect to \^f v ) 
with Ej and E v their respective energies. Similarly, the 
scalar quadrupole polarizability of the valence state \*& v ) 
is given by 

K*t,|QI*/>l a 



a°(J„) = £■ 



Ijtv 



Et — E,„ 



2 v iaii9iiJj)i 2 f241 

" 5(2 J v + 1) 2f E x -E v { ■ } 

where Q is the quadrupole operator. 

Extending these definitions, the scalar polarizability 
between two (possibly different) states \&f) and |^,) is 
given by [301 






1 



1 



Ef — Ej Ei — Ej _ 
x(* / |0 (fe) |*/)(*/|0 (fe) |* J ), 



where O^ represents the dipole operator D for k = 1 
and the quadrupole operator Q for k = 2, respectively. 
As a special case the scalar polarizabilities of a state can 
be recovered by setting i = f in the above equation. 
Apart from the static polarizability, a dynamic polariz- 
ability can also be defined. The imaginary part of the 
dynamic polarizability between two states is given by 



a°M 



E 



Ef - Ej 



Ei — Ei 



(Ef-Ej^ + uj 2 (Ei- 
x(*/|O w |*j)(*j|oW|*i), 



(2.5) 



where u is the frequency of the external electromagnetic 
field. 

From these definitions it follows that the determination 
of the polarizabilities requires the evaluation of transi- 
tion matrix elements and the excitation energies, hence 
a powerful many-body approach is necessary to evaluate 
the above quantities to high accuracy. 



B. Van der Waals coefficients 

The general expression for the van der Waals coefficients 
between two different atoms a and b in terms of their 
dynamic polarizabilities is given by [9j 



a 



ab _ (2n-2)\ 

2tt 



n-2 °° 



where V = n — I — 1 and af(ico) and a\,{\uj) are the 2 l - 
pole polarizability of atom a and 2 l -pole polarizability 
of atom b, respectively. In this article, we evaluate the 
Cq and Cs coefficients for the s-wave ground state of the 
Li atom using the simple formulas 

oc 

C 6 = ^-fdujla^itu)} 2 , (2.7) 



OC 

C* 8 = - fdw[ai(iw)a 2 (iw)], (2.8) 



obtained from Eq. (|2.6[) . The long-range part of the in- 
teraction between three ground-state atoms is not exactly 
equal to the interaction energies taken in pairs. There is 
an extra term which comes from the third-order pertur- 
bation. This correction to the van der Waals potential 
can be given as V(R) oc —v/R 3 , where [8| 



r = — I duj[ai(iuj)] , 
■k ' 
o 



(2.9) 



is called the triple-dipole constant. We have also deter- 
mined this quantity v for the Li atom and present the 
result here. 



III. METHOD OF CALCULATION 



The aim of this work is to evaluate Eq. (|2.5[) for 
both static (lo = 0) and dynamic (finite u>) polarizabili- 
ties, while avoiding the sum-over-intermediate-states ap- 
proach and at the same time treating electron-correlation 
effects rigourously. Coupled-cluster (CC) theory is one of 
the most powerful methods to incorporate the electron- 
correlation effects to all orders in the atomic wave func- 
tions. We employ here a relativistic CC theory that can 
determine the atomic wave functions accurately. 

Using Eq. (|2.5j) . we write for the dynamic polarizabil- 
ity between states |^/) and \^i) 

a*(iw) = <%|O«|*J) + <*/|0 (fc) l^>- (3.1) 



Comparing Eq. (|2.5p and Eq. (|3.1[) . we can express \^>' v ), 
where v = i, /, as 



K) = E 



Ei-Ek 



Ijtv 



(Ej - E v y 

Hi — E v 



|tfj><*/|0< fc >|*„) 



(Hj - E v y 

H — E v 
H - E v + iu 






H — E v 



(k) \^ v ) 



H - E v - 
1 



H — E„ 



H - E„ 



, (fc) 



I*.), 



(3.2) 



where H is the Dirac-Coulomb Hamiltonian. If we next 
define an effective Hamiltonian 

H eS = H - \lu 

and an effective dipole or quadrupole operator 

. (k) 



(fe ) = H - E v 
eff H- E„ + 



we can find \%f' v ) as the solution of 

(ff cff -£„)|<)=0f>„), 



(3.3) 



where \*&' v ) are the first-order perturbed wave functions 
due to the external field. 



In CC theory, the atomic wave function \^> v ) due to the 
real part of the effective Hamiltonian of a sin gle v alence 
(v) open-shell system can be expressed as [3l|, [32, l33| 



|*„)= c T {l + S v }\<S> v ), 



(3.4) 



where we define the reference state |$„) = a l\®o), with 
$o) the closed-shell Dirac-Fock (DF) state, which is 
taken as the Fermi vacuum. T and S v are the CC excita- 
tion operators for core to virtual electrons, and valence- 
core to virtual electrons, respectively. The curly bracket 
in the above expression represents the normal-ordered 
form. In our calculation, we consider all possible single 
(S) and double (D) excitations, as well as the most im- 
portant triple (T) excitations, an approximation known 
as the CCSD(T) method Hi]. To determine the ampli- 
tudes of the CC excitation operators we use 



<$ L |H c |$o) 
(**\H C S V \* V ) 



AE S L ,o, 

-<*?l#c|*„) 

+{$* \S v \$ v ){$ v \H e {l + S v }\$ v )8 K ,o 

= -($?\H C \$ V ) + ($x\S v \$ v )AE v 6 K ,o 

(3.5) 

where we have defined H c = {i?Are T } . The superscript 
L (= 1,2) represents the singly or doubly excited states 
from the closed-shell reference (DF) wave function and 
A-E is the correlation energy for the closed-shell system. 
Further, AE V is the electron affinity energy of the valence 
electron v, K (= 1,2) denotes the singly or doubly excited 
states from the single valence reference state, and the 
subscripts N and c represent the normal-ordered form 
and connected terms, respectively. Eqs. (|3.5[) are non- 
linear, and they are solved self-consistently by using a 
Jacobi iterative procedure. With the amplitudes of the 
CC excitation operators known, the zeroth-order wave 
functions can be calculated by using Eq. (|3.4[) . 



B. Determination of the first-order wave functions 

The next step is to determine the first-order wave func- 
tions. We write the wave function of a state with valence 
electron v in the presence of an external field as 



!*«> 



I*, 



i*;), 



(3.6) 



where \^ v ) is the wave function of the system in the 
absence of the external field and \^f' v ) is the first-order 
correction to |\fr„) due to the external field. In the spirit 
of the CC approach, we take the ansatz 



A. Determination of the DC wave functions 

To carry out our calculations, we will use CC cluster 
theory. As this has been described in detail in many 
other papers , we will limit ourselves to a short overview. 



|^)=e T {l + S„}|^}, 


(3.7) 


where T and S v are defined as 




f = T + T', 


(3.8) 


Oq) = Oy ~T Oy 


(3.9) 



TABLE I: The static dipole polarizability ot\ of many low-lying levels in Li [au] 



Level 


Experiments 


Other theoretical works 


This 


work 




Scalar 


Tensor 


Scalar 


Tensor 


Scalar 


Tensor 


2s 2 5i/ 2 


164(3.4)° 
164.2(1. 1) 6 




162.3% 
164 7 




162.87 




3s 2 S 1/2 


- 


4136 c 


4133% 4098^^ 
3832 d 




4107 




4s 2 S 1/2 






3.526xl0 4 % 
35040 / 




3.449 xlO 4 




2p 2 Pi/ 2 


127(3.4)' 
126.9(6) 9 




117.8° 




129.41 




3p 2 Pi/ 2 


- 




2.835xl0 4e 




2.938 xlO 4 




4p 2 Pi/ 2 


- 




2.734xl0 5e 




2.635 xlO 5 




2p 2 P 3 /2 


127.2(7) 9 


1.64(4) 9 


117.8 e 


3.874 e 


123.09 


5.95 


3p 2 P 3 /2 


- 


- 


2.835xl0 4e 


-2173 e 


2.929 xlO 4 


-2078 


4p 2 P 3 /2 


- 


- 


2.735xl0 5e 


-2.074xl0 4e 


2.634 xlO 5 


-1.473 xlO 4 


3d 2 Z> 3/2 


-15130(40) h 


1.643(6) xl0 4h 


-1.504xl0 4e 


1.147xl0 4e 


-1.953 xlO 4 


1.412 xlO 4 


4d 2 D 3/2 


- 


- 


3.093xl0 6e 


-5.355xl0 5e 


3.834 xlO 6 


-6.650 xlO 5 


3d 2 L» 5/2 


-15130(40)'' 


- 


-1.510xl0 4e 


1.645xl0 4e 


-2.008 xlO 4 


2.139 xlO 4 


4d 2 £> 5 /2 


- 


- 


3.103xl0 6e 


-7.678 xl0 5e 


3.843 xlO 6 


-9.496 xlO 5 



a Molof et al. (1974) [35j, " Miffre et al. (2006) [36j, c Themelis et al. (1995) [22J], 
Ashby et al. (2003) 0, / Magnier et al. (2002) [2j|, 9 Windholz et aMl992) [l| ( 6 Li 

i Hunter et al. (1991) [lj. 



d Merawa et al. (1998) [23J 
values), h Ashby et al. (2003) [39 



Here T" and S' v are the corrections to the T and S v oper- 

(k) 

ators in the presence of the operator 0^ ff , respectively. 
Substituting Eqs. flU]) and flU) in Eq. <J3TTI) , we find 

|tf„) - e T [l + S v + r{l + S v } + S' v }\<P v ), 

(3.10) 

where only the terms linear in T" and S' v exist, since Eq. 
contains just one 0^ ff operator. By comparing Eqs. 



(IQl . (EU), and (l3~T0l) . we get 

K) = e T [T'{l + ^} + 5;,]|^). (3.11) 

We evaluate these perturbed CC operator amplitudes us- 
ing the following equations (c/. Eqs. (|3.5[) ) : 

($ L | [H* + u 2 ]T'\$ ) = ($ L \ {()We T \ |$o), 



C. Evaluation of a using the RCC approach 

The expression for the polarizabilities using our CC ap- 
proach can now be obtained by substituting Eqs. (|3.4[) 
and (|3.11[) in Eq. (|3. 1 [) . In this way we get (we also 
normalize) 



afe(iw) 






(<*,|{1 + S}}~m[T'{l + Si} + agi*,>+ 

($/|[Sjf + {1 + S}}T't]OW{l + S«}|* 4 >) 

(3.13) 



where 



($f|[(tf c - AE v ) 2 +lu 2 ]S' v \<I> v ) 

- <$f|^ c + ^^r( 1 ){i + 5(°)}|<i>„) 

(3.12) 

where the meaning of L and K was explained above. 
The first-order wave functions are determined using Eq. 
p. lip after obtaining the perturbed CC amplitudes. 



N v = <$„|{1 + SlWo{l + S v }\$ v ), 

with v — i, f, and we have defined 0( fc ) — e T 0^e T and 
JV„ = e Tt e T . 

We first evaluate, by using the generalized Wick's the- 
orem, the intermediate terms Ow and A/"o in the above 
expressions as effective one-body, two-body, and so on, 
terms. Next we sandwich the open-shell valence-core 
electron excitation operators to evaluate the exact ex- 
pression. 



IV. RESULTS AND DISCUSSIONS 

We have used partly numerical and partly analytical or- 
bitals to generate the complete basis sets. The numerical 
orbitals were obtained using GRASP 41| , and the analyt- 
ical orbitals were obtained using Gaussian-type orbitals 
(GTO's) [13]. In total, we have taken up to the 30s, 30p, 
25d, 25 /, and 20g orbitals to calculate the DF wave func- 
tion. Out of these, we have generated the first 4, 3, 2, 2, 
and 2 orbitals from the s, p, d, /, and g symmetries, re- 
spectively, using GRASP. The remaining continuum or- 
bitals were obtained analytically from GTO's, using as 
parameters a — 0.00525 and (3 — 2.73. After this, the 
final orbitals were orthogonalized using Schmidt's proce- 
dure [H. 



f T ;u) 





-oo 



(ii) 



FIG. 1: Extra correlation diagrams which appear in the calcu- 
lation of the polarizabilities using our novel approach. These 
diagrams do not appear when the CC wave functions are used 
in the sum-over-the-states method. 

We present the static dipole and quadrupole polariz- 
abilities of several important low-lying states of Li in Ta- 
ble U and Table UH respectively. In these Tables, we have 
also listed other theoretical results and the most recent 
experimental results, where available. For the ground 
state, a number of theoretical dipole polarizability results 
are available, for the excited states, however, few calcula- 
tions have been carried out. All other theoretical results 
except one are based on non-relativistic theory. Some 
of these calculations are also performed using molecular 
codes, at the cost of atomic symmetries [25]. The one 
available relativistic calculation on the excited states is 
carried out using a rather approximate method to in- 
clude the correlation effects due to the Coulomb inter- 
action 37]. Our calculation uses a relativistic approach 
which considers correlation effects to all orders in the 
form of CC amplitudes. Table HT1 shows the result for the 
static quadrupole polarizabilities. No experimental data 
is available for comparison, and the available theoretical 
results for the 2S l 1 / 2 level are not very consistent. 

Although our method is theoretically superior to the 
previously employed methods to determine both dipole 
and quadrupole polarizabilities, it seems that some of 
the earlier results are in better agreement with the ex- 
perimental results than ours. This may be due to the 
fact that experimental energies are used in some of these 
calculations in contrast to our method which is fully ab 
initio. This means that in our calculation there may 



L 2 P 

* 3/2 


oo 


v N V 


\ 


*2p 

3/2 


^o 



FIG. 2: The correlation diagram that causes a large discrep- 
ancy between the calculated and the experimental results of 
the tensor polarizability of the 2p 3 / 2 state. 



TABLE II: The static quadrupole polarizability a® of many 
important states in Li [au]. 



Level 


Other theoretical works 


This work 


2s 2 S 1/2 


1423", 
1403 e 


1424", 1430", 1423 
, 1393 / , 1424(4) 9 , 


.266(5) d , 
1424.4 /l 


1420 


3s 2 S 1/2 




3.5642xl0 5h 




3.475 xlO 5 


4s 2 S 1/2 




1.1587xl0 7h 




1.113xl0 7 


2p 2 Pi/ 2 




- 




7.804 xlO 4 


3p 2 Pi/ 2 




- 




1.033 x 10 7 


4p 2 Pi/ 2 




- 




3.301 x 10 9 



Spelsberg et al. (1993UM, b Marinescu et al. (1994) [11], 
Merawa et al. (1994) iff" Yan et al. (1996) [|, e ' f Patil 
and Tang (1997,1999) g(| S3, 9 Snow et al. (2005) [3, 
h Zhang et al. (2007) [13]. 



be strong cancelations with neglected higher-order exci- 
tations in the correlation effects. We note that in our 
approach we implicitly take into account certain corre- 
lation effects that cannot be accounted for in the usual 
sum-over-states approach that is used in so many of the 
earlier calculations. These diagrams, which are shown 
diagrammatically in Fig. [TJ are part of the RPA. 

As Table U shows, our value for the tensor polarizabil- 
ity of the 2p 2 P 3 / 2 level is larger than the experimental 
result. In our investigation we found that this large value 
is due to the unusual behavior of the correlation effects 
produced by the diagram shown in Fig. O Leaving out 



TABLE III: The off-diagonal scalar polarizability in Li [au] 





DF 




CCSD(T) 


Dipole 








2s -3s 


-27.18 




-20.41 


2s -4s 


-202.9 




-164.2 


3s -4s 


-105.8 




6.292 


Quadrupole 








2s -3s 


2.495 x 


It) 4 


2.219xl0 4 


2s -4s 


1.245 x 


10 5 


1.134 xlO 5 


3s -4s 


9.281 x 


10" 


6.647 xlO 5 




o , ccsd(7)-(1'df lau] (dipole) 



a Q (i u) [auj 


1400 




1200 




1000 




800 




600 




400 




200 














FIG. 3: The imaginary parts of the dipole (i) and quadrupole 
(ii) polarizabilities of the ground state of Li as a function of 
the angular frequency u). 



TABLE IV: Cs and C& values for the ground states of Li-Li 
faul. 



C 6 (xl0 3 ) C 8 (xl0 5 ) 



This work 

Dirac-Fock 

CCSD(T) 



1.473 0.8891 

1.396(6) 0.8360 



1.388 



Other theoretical works 

Marinescu et al. (1994) [44] 

Spelsberg et al. (1996) [16] 

Yan et al. (1996) [8] 

Patil and Tang (1999) [47] 1.360 

Porsev and Derevianko (2003) [49] - 

Mitroy and Bromley (2003) [18] 1.3946 



0.8324 



1.39322 0.834258(42) 
0.8100 
0.834(4) 
0.83515 



this diagrams yields a value for the tensor polarizability 
of the 2p 2 P 3 / 2 level of ~ 1.6, which agrees nicely with 
the experiment. For the completeness of the theory this 
effect cannot be left out. We expect that this effect will 
cancel with the neglected higher-order excitations. 

In Table Hill we present scalar dipole and quadrupole 
polarizabilities among different s-states of Li which are 
also important in the determination of the van der Waals 
coefficients of the excited states for ultra-cold atom ex- 
periments. Our method can also be employed to deter- 
mine these quantities in the heavy alkali atoms like Cs 
and Fr that are important candidates for the study of 
atomic parity nonconservation [27]]. To our knowledge, 
no other results are available to compare with these re- 
sults. As the Table shows, the scalar dipole polarizability 
between the 2s and 3s states in Li is of opposite sign to 




o.o 



0.1 



0.2 



0.3 



0.4 



d [auj 



(0 



accsDm-ffDF lau] (quadrupole) 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 
(ii) 



l> [auj 



FIG. 4: The difference between the DF and CCSD(T) results 
for the imaginary parts of the dipole (i) and quadrupole (ii) 
polarizabilities of the ground state of Li as a function of the 
angular frequency u. 



TABLE V: The triple-dipole constant v for Li-Li-Li (xlO 4 
faul. 



V (Li-Li-Li) 



This work 

Dirac-Fock 

CCSD(T) 

Other theoretical works 

Yan et al. (1996) [8] 

Mitroy and Bromley (2003) [18] 



18.576 
16.934 



17.0595(6) 
17.087 



the other alkali atoms [271 ]. 

The main goal of this work is to illustrate how to 
evaluate the van der Waals coefficients using the present 
method. Figure [3] shows the imaginary parts of the dipole 
and quadrupole polarizabilities (in atomic units) of the 
ground state of Li as functions of angular frequency, to. 
As the Figures show, these quantities fall off exponen- 
tially for higher values of ui. To illustrate the effect of 
electron correlation as a function of frequency, we have 
plotted the difference between the CCSD(T) and the DF 
results in Fig. 2] This Figure suggests that the corre- 
lation effects vanish for higher frequencies. Using the 
imaginary parts of the dipole and quadrupole polariz- 
abilities in Eqs. ([2~7) . ([2~8]) and ([2~9]) . we evaluated the 
Cq, C's, and v coefficients, respectively, using a numerical 
integration method. 

In Table HVl we present our C@ and Cs coefficients and 
compare them with the other available results. Although 



our value for the static polarizability of the ground state 
of Li is slightly smaller than the results presented by 
others, our C$ and Cg values are in good agreement with 
the other results. We present the coefficient v of the 
third-order correction to the long-range potential in Ta- 
ble [V] which matches well with the other available semi- 
empirical results. 



of the polarizabilities which we used to evaluate the 
van der Waals coefficients of the Li atom. By us- 
ing this novel technique, we were able to consider the 
electron-correlation effects rigorously because the tech- 
nique is fully relativistic and it avoids the sum-over-states 
method. 
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